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ABSTRACT 

High-resolution  measurements  obtained  from  NOAA  “best”  atmospheric  turbulence  (BAT)  probes 
mounted  on  an  EGRETT  high-altitude  research  aircraft  were  used  to  characterize  turbulence  in  the  upper 
troposphere  and  lower  stratosphere  at  scales  from  2  m  to  20  km,  focusing  on  three-dimensional  behavior  in 
the  sub-kilometer-scale  range.  Data  were  analyzed  for  129  separate  level  flight  segments  representing  41  h  of 
flight  time  and  12  600  km  of  wind-relative  flight  distances.  The  majority  of  flights  occurred  near  the  tropopause 
layer  of  the  winter  subtropical  jet  stream  in  the  Southern  Hemisphere.  Second-order  structure  functions  for 
velocity  and  temperature  were  analyzed  for  the  separate  level-flight  segments,  individually  and  in  various 
ensembles.  A  3D  scaling  range  was  observed  at  scales  less  than  about  100  m,  with  power-law  exponents  for  the 
structure  functions  of  the  velocity  component  in  the  flight  direction  varying  mostly  between  0.4  and  0.75  for  the 
separate  flight  segments,  but  close  to  %  for  the  ensemble-averaged  curves  for  all  levels  and  for  various  sub¬ 
ensembles.  Structure  functions  in  the  3D  scaling  range  were  decoupled  from  those  at  scales  greater  than  10  km, 
with  the  large-scale  structure  functions  showing  less  variation  than  those  at  smaller  scales.  Weakly  anisotropic 
behavior  was  observed  in  the  3D  range,  with  structure  parameters  for  the  lateral  and  vertical  velocities  on  the 
same  order  as  those  in  the  flight  direction  but  deviating  from  the  expected  isotropic  value.  Anisotropy  was 
correlated  with  turbulence  intensity,  with  greater  anisotropy  associated  with  weaker  turbulence. 


1.  Introduction 

Turbulence  in  the  upper  troposphere  and  lower  strato¬ 
sphere  (UTLS)  can  adversely  impact  the  performance  of 
a  range  of  aerospace  systems.  Optical  turbulence  (OpT) 
from  fluctuations  in  temperature  and  humidity  can  affect 
electromagnetic  propagation,  disrupting  communications, 
radar,  high-energy  laser  systems,  and  space  imaging. 
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Similarly,  clear-air  turbulence  (CAT)  can  lead  to  aircraft 
upset,  a  problem  particularly  acute  for  stratospheric  ve¬ 
hicles  that  often  must  fly  within  a  narrow  flight  envelope. 
Thus,  reliable  prediction  is  critical  for  turbulence  avoid¬ 
ance  and  mitigation  in  these  operational  situations.  Un¬ 
fortunately,  prediction  is  difficult  because  of  the  nature 
of  UTLS  turbulence,  which  is  concentrated  in  very  thin 
layers  with  thicknesses  of  hundreds  of  meters  to  a  few 
kilometers,  a  scale  smaller  than  the  typical  resolution  of 
mesoscale  weather  forecasting  models.  As  such,  pre¬ 
diction  and  forecasting  rely  on  parameterized  models, 
based  on  statistical  correlations  and  phenomenological 
relations  between  mesoscale  parameters  and  turbulence 
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quantities,  in  order  to  identify  regions  that  have  the  po¬ 
tential  for  high  levels  of  clear  air  or  optical  turbulence. 
These  approaches  are  only  as  good  as  the  fundamental 
understanding  of  UTLS  turbulence  structure  that  forms 
the  underlying  foundation  for  the  models;  such  insight 
must  be  gleaned  from  experimental  data,  detailed  mi¬ 
croscale  simulations,  highly  nested  mesoscale  model  runs, 
or  a  combination  of  all  three. 

Two  seminal  experiments  have  formed  the  primary 
database  for  studying  UTLS  turbulence,  both  involving 
instrumentation  of  commercial  transport  aircraft.  The 
Global  Atmospheric  Sampling  Project  (GASP)  utilized 
specially  selected  Boeing  747  in  routine  commercial  air¬ 
line  service  from  1975  to  1979,  representing  6900  flights, 
instrumented  to  measure  horizontal  wind  velocity,  tem¬ 
perature,  and  trace  constituents  (Perkins  1976;  Nastrom 
and  Jasperson  1983).  A  summary  of  the  GASP  dataset 
can  be  found  in  Jasperson  et  al.  (1990)  and  Nastrom  and 
Gage  (1985).  The  Measurement  of  Ozone  and  Water 
Vapor  by  Airbus  In-Service  Aircraft  (MOZAIC)  project 
utilized  five  specially  instrumented  Airbus  aircraft  in 
normal  commercial  operation  from  1994  to  1997,  repre¬ 
senting  over  7500  flights  and  54  000  flight  hours  (Marenco 
et  al.  1998).  Although  the  major  emphasis  of  the  program 
was  water  vapor  and  ozone  measurements,  meteorolog¬ 
ical  data  included  horizontal  wind  velocity  and  temper¬ 
ature,  with  data  sampling  at  4-s  intervals  on  all  flights. 

The  wealth  of  data  from  these  two  campaigns,  and  the 
reinterpretation  of  the  data  by  Lovejoy  et  al.  (2009)  and 
Lilley  et  al.  (2008),  has  furthered  the  understanding  of 
UTLS  turbulence;  however,  spatial  resolution  was  lim¬ 
ited  due  to  sampling  rates  and  the  speed  of  the  aircraft, 
restricting  investigation  of  turbulence  structure  and 
phenomena  to  larger  scales — from  a  few  kilometers  up 
to  thousands  of  kilometers.  Unfortunately,  this  range  of 
scales  does  not  encompass  the  key  features  of  CAT  or 
OpT,  both  of  which  are  characterized  by  scales  in  the  1  m 
to  1  km  range  and  are  associated  with  three-dimensional 
phenomena  that  cannot  be  investigated  using  only  hori¬ 
zontal  velocities  (e.g.,  the  vertical  velocity,  which  is  a 
critical  contributor  to  aircraft  response  to  CAT). 

Although  aircraft  data  at  scales  down  to  about  80  m 
are  available  from  other  aircraft  campaigns,  these  are  as¬ 
sociated  with  limited  numbers  of  reported  flight  segments 
(e.g.,  Kennedy  and  Shapiro  1980;  Chan  et  al.  1998). 

One  solution  to  the  scale  problem  has  been  to  extrap¬ 
olate  large-scale  experimental  or  mesoscale  model  data 
to  microscale  behavior.  Lindborg  (1999,  hereafter  L99) 
and  Cho  and  Lindborg  (2001,  hereafter  CL01)  estimated 
dissipation  from  MOZAIC  data  by  fitting  structure  func¬ 
tions  and  assuming  two-dimensional  turbulence  behav¬ 
ior.  Frehlich  and  Sharman  (2005,  hereafter  FS05)  used 
the  curve  fits  from  L99  to  extrapolate  mesoscale  model 


results  to  smaller  scales  for  estimating  dissipation,  struc¬ 
ture  functions,  and  other  characteristics  of  microscale  tur¬ 
bulence.  However,  these  methods  have  not  been  verified 
because  of  the  dearth  of  available  data  at  small  scales,  as 
pointed  out  in  FS05.  This  paper  will  show  that  such  ex¬ 
trapolations  are  not  universally  valid. 

The  work  presented  in  this  paper  aims  to  fill  the  gap  in 
published  UTLS  experimental  data  for  turbulence  of 
scales  in  the  range  from  approximately  1  m  to  1  km.  The 
work  focuses  on  analysis  of  structure  functions  calcu¬ 
lated  using  high-resolution  turbulence  measurements  of 
temperature  and  all  three  components  of  velocity  as  part 
of  a  multiyear  effort  (1998-2006)  to  study  CAT  and 
OpT.  The  overall  objective  is  to  probe  the  small-scale 
structure  of  UTLS  turbulence  in  the  hopes  of  further¬ 
ing  our  understanding  of  atmospheric  phenomenon  in 
general  and  of  providing  data  for  building  and  testing 
more  effective  parameterized  models  for  prediction  and 
forecasting. 

2.  Measurements 

Aircraft  measurements  were  acquired  with  the  Grob 
520T  EGRETT  aircraft,  the  sole  twin-seat  version  of  the 
four  Grob  520  aircraft  produced,  which  was  operated  by 
Airborne  Research  of  Australia  of  Flinders  University. 
This  high-altitude,  turboprop  research  aircraft  is  capable 
of  operation  at  altitudes  up  to  15  km  at  air  speeds  of  ap¬ 
proximately  100  m  s  1  with  an  endurance  of  8  h.  The  air¬ 
craft  can  accommodate  up  to  three  National  Oceanic  and 
Atmospheric  Administration  (NOAA)-built  “best”  atmo¬ 
spheric  turbulence  (BAT)  probes  (Crawford  and  Dobosy 
1992,  1997),  one  located  under  each  wing  and  one  at  the 
top  of  the  tail.  The  structure  functions  obtained  from  the 
two  underwing  probes  were,  in  general,  nearly  identi¬ 
cal,  so  the  left  BAT  probe  was  used  for  analysis,  except  in 
cases  when  its  signal  was  noisier  than  that  from  the  right 
probe  or  for  the  occasional  flights  when  data  were  missing 
or  out  of  range.  The  tail  probe  data  were  not  used  in  this 
study. 

The  BAT  probe  has  a  13-cm-diameter  hemisphere  with 
nine  pressure  ports  for  measurement  of  all  three  com¬ 
ponents  of  wind  velocity  and  a  microbead  thermistor 
located  inside  the  central  dynamic  pressure  port  for  tem¬ 
perature  measurement.  Beginning  in  2002,  the  NOAA 
Fast  Ultrasensitive  Temperature  (FUST)  probe  was  also 
mounted  externally  on  the  top  of  the  BAT  probe  body. 
This  bare-wire  microthermocouple  sensor  featured  in¬ 
creased  signal-to-noise  compared  to  the  thermistor,  so  its 
output  was  used  when  available.  For  three  of  the  cam¬ 
paigns  (described  below),  the  right-wing  BAT  probe  was 
replaced  with  a  standard  Rosemont  five-hole  probe  with 
a  Rosemont  PT50  probe  for  temperature. 
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Wind  velocities  were  obtained  by  subtracting  pressure- 
sphere-derived,  wind-relative  aircraft  velocity  measure¬ 
ments  from  the  ground-fixed  aircraft  velocities  and 
orientations  (pitch,  yaw,  and  roll)  obtained  from  the  In¬ 
ertial  Navigation  System  (INS)  and  GPS  measurements. 
With  this  approach,  vertical  wind  velocity  is  most  prone 
to  error,  since  it  is  found  by  subtracting  two  numbers  that 
are  nearly  equal  to  obtain  a  small  number,  and  because 
it  is  more  sensitive  to  aircraft  angle  of  attack  than  the 
horizontal  wind  velocities.  Garman  et  al.  (2006)  found  a 
lcr  precision  of  0.03  m  s  1  for  the  vertical  velocity  mea¬ 
sured  with  the  BAT  probe  after  being  calibrated  in  a 
wind  tunnel.  However,  without  such  calibration,  the  over¬ 
all  error  in  the  wind  components  is  closer  to  0.25  m  s  1 
(Crawford  and  Dobosy  1992).  For  comparison,  Nastrom 
et  al.  (1984)  report  an  overall  accuracy  of  0.5  m  s  1  for 
the  GASP  measurements  and  in  L99,  a  similar  level  of 
accuracy  was  reported  for  the  MOZAIC  measurements. 

For  nearly  all  flights  prior  to  2006,  velocity  and  tem¬ 
perature  data  were  sampled  at  50,  55.6,  or  58.8  FIz,  pro¬ 
viding  horizontal  spatial  resolution  of  approximately  1.4 
to  2  m.  For  2006,  the  sampling  rate  was  10  Hz  and  for 
November  2002,  it  was  25  Hz,  reducing  spatial  resolution 
to  approximately  10  m  and  4  m,  respectively.  Frequency 
response  of  the  velocity  measurements  was  generally 
near  or  above  the  sampling  frequency.  On  the  other  hand, 
the  time  constant  of  the  BAT  thermistor  probe  resulted 
in  a  thermal  lag  drop-off  at  scales  near  20  to  30  m.  Both 
the  FUST  and  the  Rosemont  temperature  sensors  had 
slightly  faster  response,  but  still  exhibited  drop-off  in  the 
10-  to  20-m  scale  range. 

A  typical  flight  featured  several  level  flight  segments 
at  altitudes  from  7  up  to  14  km,  covering  wind-relative 
distances  from  45  to  350  km.  The  EGRETT  was  flown 
on  autopilot,  set  to  maintain  constant  pressure  altitude, 
with  altitude  variations  less  than  ±25  m  for  all  but  five 
of  the  129  flights.  As  a  general  rule,  the  segments  were 
flown  upwind  or  downwind,  with  over  75%  of  the  dis¬ 
tances  flown  with  angles  between  wind  and  flight  path  of 
less  than  30°.  Approximately  5%  of  the  flight  distance 
involved  significant  crosswind  components,  with  angles 
greater  than  70°.  Data  were  also  available  for  climb  and 
descent  between  level  segments  for  some  flights,  and 
these  were  used  to  estimate  vertical  gradients. 

Lovejoy  et  al.  (2004,  hereafter  LST04)  have  shown  that 
spectra  and  structure  functions  derived  from  airborne 
turbulence  measurements  can  be  strongly  influenced  by 
aircraft  trajectories,  an  effect  due  to  anisotropic  scaling  in 
the  horizontal  and  vertical  directions,  combined  with 
vertical  excursions  of  the  aircraft  in  response  to  the  tur¬ 
bulence  motions.  This  effect  is  seen  at  large  scales,  where 
aircraft  inertia  no  longer  attenuates  altitude  fluctuations. 
Following  the  approach  in  LST04,  the  extent  of  the  aircraft 


Fig.  1.  First-order  structure  function  of  aircraft  altitude  fluctu¬ 
ations  for  1  Sep  1998  at  13.1  km.  Dashed  line  is  Dz  =  sr,  where  s  is 
the  slope  of  the  aircraft  trajectory. 


inertia  range  was  found  from  the  first-order  structure 
function  of  aircraft  altitude  fluctuations: 

Z)z(r)  =  (|zAc(x)“zAc(*  +  ,')l)-  (x) 

Note  that  Dz  follows  a  linear  scaling  with  r  in  the  aircraft- 
inertia  range,  estimated  to  be  at  scales  less  than  1  km,  as 
seen  in  Fig.  1  for  one  of  the  levels.  The  1-km  upper  limit 
is  slightly  lower  than  the  ER-2  value  of  3  km  (LST04) 
because  of  the  lower  wing  loading  of  the  EGRETT 
(25  lb  ft~2  compared  to  40  lb  ft  2).  The  Dz  plots  for  all 
other  cases  showed  aircraft -trajectory  slopes  on  the  order 
of  1  to  4  m  km  ,  similar  to  values  seen  in  ER-2  data 
(LST04)  and  estimates  for  the  MOZAIC  data  (Lilley 
et  al.  2008).  Based  on  these  results,  quantitative  analysis 
of  structure  functions  will  be  confined  to  sub-kilometer 
scales,  and  interpretation  of  larger-scale  behavior  will  be 
limited  to  qualitative  observations.  Additional  discus¬ 
sion  of  the  effects  of  aircraft  trajectories  on  measured 
structure  functions  can  be  found  in  section  5a. 

The  129  level  flight  segments  selected  for  analysis  are 
summarized  in  Table  1.  They  represent  41  h  of  flight  time 
and  12  600  km  of  wind-relative  flight  distances,  much  less 
than  the  MOZAIC  and  GASP  campaigns,  but  significant 
given  the  higher  data  rates  and  the  flight  planning  re¬ 
quired  for  targeting  strong  turbulence  periods. 

To  maximize  the  likelihood  of  finding  turbulence,  the 
measurement  flights  were  focused  on  the  winter  sub¬ 
tropical  jet  stream,  primarily  in  and  around  Adelaide, 
Southern  Australia,  which  was  the  base  of  operation  for 
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Table  1.  Summary  of  level  flight  segments  (SA,  South  Aus¬ 
tralia;  JA,  Japan;  WA,  Wales;  DA,  Darwin,  Australia;/),  sampling 
frequency;  D,  total  distance  of  level  segments  in  wind-relative 
coordinates). 


Dates 

Loc. 

No. 

segments 

Alt  (km) 

fs  (Hz) 

D  (km) 

1998:  8/25-9/03 

SA 

25 

6.8-14.2 

55.6 

1900 

1999:  2/03-2/21 

JA 

19 

5.8-12.4 

58.8 

2026 

1999:  8/6 

SA 

7 

9.0-12.2 

58.8 

710 

2000:  5/10-6/6 

WA 

20 

8.1-13.8 

55.6 

2112 

2001:  8/17-8/24 

SA 

15 

8.5-10.8 

55.6 

1335 

2002:  8/28-9/16 

SA 

37 

7.1-14.0 

50 

4101 

2002:  11/22 

DA 

5 

9.6-14.4 

25 

442 

2006:  8/6 

SA 

1 

8.5 

10 

50 

the  EGRETT.  Japan  was  the  location  for  one  of  the  two 
1999  campaigns,  with  the  goal  of  probing  the  strong 
Northern  Hemisphere  winter  subtropical  jet  stream  in 
that  region.  The  measurements  from  the  Adelaide  and 
Japan  flights  represent  80%  of  the  total  flight  distances. 
The  remaining  flights,  in  May  and  June  of  2000  and  in 
November  2002,  were  associated  with  two  projects  that 
utilized  the  EGRETT  for  primary  purposes  other  than 
CAT  or  OpT  in  and  near  jet  streams.  The  AberEgrett 
Experiment  (Whiteway  et  al.  2003a,  2004)  involved  a 
consortium  of  UK  researchers,  and  was  conducted  over 
southwest  England  in  summer  of  2000  to  study  gravity 
waves  induced  by  the  Welsh  mountains  and  their  effects 
on  the  tropospheric-stratospheric  exchange  of  trace  gases. 
Similarly,  the  EMERALD  2  (Whiteway  et  al.  2003b) 
campaign  was  conducted  over  Darwin,  Australia,  during 
the  winter  of  2002  to  study  cirrus  clouds.  These  two  cam¬ 
paigns  employed  the  Rosemont  probe  under  the  right 
wing. 

Histograms  of  several  flight  parameters,  weighted  based 
on  the  length  of  each  segment  in  wind-relative  distance, 
are  shown  in  Fig.  2.  The  mean  altitude  was  10.4  km  and 
the  mean  wind  speed  was  44.5  m  s-1,  both  exhibiting 
Gaussian-like  distributions.  The  wind  direction  statistics 
reflect  prevailing  westerly  jet  streams,  with  a  mean  direc¬ 
tion  of  254°  and  with  84%  of  the  flight  distances  within 
±30°  of  wind  from  due  west.  The  flight  segment  lengths 
exhibit  a  skewed  distribution,  with  more  than  75%  of  the 
flight  distance  flown  in  segments  less  than  the  mean  of 
123  km.  The  tropopause,  defined  here  as  the  lowest  al¬ 
titude  for  which  the  lapse  rate  was  less  than  2°C  knU1 
[the  World  Meteorological  Organization  (WMO)  defi¬ 
nition],  was  above  11  km  for  all  flight  days  and  above 
13  km  for  85%  of  the  flights.  These  high  tropopause 
levels,  found  using  nearby  atmospheric  soundings  (http:// 
weather.uwyo.edu/upperair/sounding.html),  are  not  un¬ 
typical  for  midwinter  in  southern  Australia,  resulting  in 
95%  of  the  flight  distances  below  the  tropopause. 
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Fig.  2.  Histograms  showing  distributions  of  flight  parameters. 
Dashed  lines  mark  mean  values. 


3.  Structure  functions 

Structure  functions  were  calculated  from  time  series 
of  velocity  and  temperature: 

^LlW  =  <[“ L +  WV,TAs)  -  UM) 

£>NN W  =  <[ UN 0  +  r/VTAs)  “  UN(0f) 

^wwO)  =  (K*  +  WFtas)  “  w0)]2> 

D^ir)  =  ([T(t  +  rlVTAS)  -  T(t)}2) ,  (2) 

where  uL  is  the  component  of  the  wind  along  the  di¬ 
rection  of  the  separation  distance  vector  r,  the  direction 
of  the  flight  path  in  this  case;  uN  is  the  component  nor¬ 
mal  to  r;  w  is  the  vertical  velocity;  T  is  temperature;  Vtas 
is  the  true  airspeed  of  the  aircraft;  and  the  angle  brackets 
denote  ensemble  averaging. 

Antonia  and  Smalley  (2001)  showed  that  structure 
functions  generally  exhibit  smaller  regions  of  power-law 
behavior  than  spectra  and  Nichols-Pagel  et  al.  (2008)  in¬ 
dicated  that  multitaper  spectra  and  wavelet  variances 
were  preferred  over  structure  functions.  Despite  this, 
structure  functions  provide  some  advantages  over  power 
spectra  (L99)  and  thus  are  adopted  for  this  study.  Struc¬ 
ture  functions  are  often  harder  to  interpret  than  spectra, 
lacking  the  direct  physical  interpretation  of  power  spec¬ 
tra.  However,  structure  functions  are  similar,  though  not 
identical,  to  cumulative  spectra,  and  thus  they  can  be 
interpreted  in  a  similar  manner — i.e.,  Z)LL(r)  is  a  relative 
measure  of  the  energy  contained  in  (uj)  for  scales  <r. 
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Equations  (2)  are  based  on  a  separation  distance  de¬ 
fined  in  wind  relative  coordinates,  r  =  VtasAi,  where  At 
is  the  time  difference  between  data  points  in  the  measured 
time  series.  The  use  of  wind-relative  distances,  which  is 
consistent  with  the  frozen  flow  hypothesis  of  G.  I.  Taylor, 
differs  from  the  approach  of  CL01,  for  which  latitude 
and  longitude  were  used  to  find  r  in  an  earth-fixed  co¬ 
ordinate  frame.  However,  the  earth-fixed  approach  can 
lead  to  errors  with  slower-flying  aircraft,  best  illustrated 
by  considering  a  case  for  which  the  aircraft  is  flying  into 
a  strong  head  wind  with  wind  speed  equal  to  true  air¬ 
speed,  such  that  ground  speed  is  0.  (A  similar  situation 
occurred  in  1999  during  an  EGRETT  overwater  flight 
south  of  Japan  where  the  jet  was  strong  enough  to  cause 
a  negative  ground  speed.)  For  this  zero-ground  speed 
example,  the  aircraft  is  stationary  in  an  earth-fixed  frame, 
so  the  fixed-earth  approach  would  erroneously  lead  to 
r  =  0  for  any  two  data  points,  regardless  of  the  time  dif¬ 
ference.  With  either  approach,  nonstationary  effects  will 
lead  to  errors  at  large  separation  distances. 

Structure  parameters  were  calculated,  assuming  a 
structure-function  scale  range  with  an  r2'  dependency: 


where  X  denotes  either  L,  N,  W,  or  T. 

Isotropic  turbulence  theory  predicts  the  following 
relations: 


fC2\m 

e=(-J  ’ 

(4) 

V  «  J 

c2  =  C2  =  -C2 

(5) 

where  s  is  the  dissipation  rate  and  a  is  a  constant,  with 
atmospheric  studies  indicating  a  value  of  a  2  (Wyngaard 
and  Cote  1971). 

Average  structure  functions  were  found  for  each  flight 
segment  for  51  values  of  r,  evenly  spaced  in  logarithmic 
terms  from  2  m  to  31  km.  The  actual  separation  dis¬ 
tances  rA  for  each  segment  varied  because  of  differences 
in  true  airspeed  and  sampling  rates.  As  such,  values  at 
the  51  specified  distances  were  found  by  interpolation  of 
the  two  surrounding  values  of  rA  using  a  power-law  fit  rn. 

Scaling  exponents  were  found  for  ranges  of  r  at  sub¬ 
kilometer  scales  that  displayed  linear  behavior  on  log- 
log  plots.  This  is  illustrated  in  Fig.  3  for  25  August  1998, 
showing  two  constant  slope  regions  demarcated  by  the 
vertical  lines.  For  each  flight  segment,  the  linear-slope 
region  with  the  exponent  closest  to  %  for  l)IA  was  used 
for  calculating  the  structure  parameters  by  force-fitting 
an  i'22'  curve.  For  example,  Cj  is  the  value  that  mini¬ 
mizes  the  mean  square  error  between  C2r2/3  and  Di  i  (r) 


Fig.  3.  Structure  functions  for  25  Aug  1998, 12.3  km;  X  refers  to 
parameter  plotted — Z)LL  (m2  s~2;  solid  line);  £>NN  (m2  s~2:  dashed 
line);  D Ww  (m2  s-2;  solid  line  with  plus  symbols);  and  Z)tt  (K2; 
dotted-dashed  line).  Vertical  dashed  lines  mark  the  beginning  of 
constant  slope  regions  and  vertical  dotted  lines  mark  the  end. 

in  the  linear  slope  region.  This  is  illustrated  graphically 
in  Fig.  3  by  the  dotted  lines  that  extend  from  each  of 
the  structure  function  curves  to  the  vertical  axis,  with  the 
intercept  at  r  =  1  marking  the  structure  parameter.  The 
resulting  exponents  for  this  case  were  very  close  to  Vv. 
0.66,  0.61,  0.69,  and  0.65  for  Z)LL,  DNN,  Dww,  and  DTT, 
respectively.  This  region  will  be  referred  to  as  the  3D 
scaling  range.  (The  adjective  “3D”  is  used  to  describe 
a  range  of  scales  for  which  the  vertical  velocity  structure 
functions  are  on  the  same  order  as  those  of  the  hori¬ 
zontal  velocities.) 

Note  that  the  /)ww  curve  in  Fig.  3  exhibits  a  leveling 
off  at  scales  above  1  km.  This  plateau  shape  is  expected 
for  separation  distances  much  greater  than  the  largest 
scale  of  the  energy-containing  eddies,  with  the  plateau 
value  equal  to  twice  the  variance  of  the  signal  (Monin 
and  Yaglom  1975).  For  measurements  that  span  a  region 
much  larger  than  the  integral  scale  of  the  flow,  structure 
functions  should  end  in  a  plateau  (referred  to  here  as  a 
terminal  plateau).  Structure  functions  that  end  with  a 
nonzero  slope  indicate  that  sample  regions  are  smaller 
than  the  largest  scales;  the  DNN  and  /)TT  curves  in  Fig.  3 
are  examples  of  the  latter.  The  shape  of  the  Dww  curve 
in  Fig.  3,  including  the  hump  preceding  the  plateau,  can 
be  generated  by  low-pass  filtering  data  that  does  not 
exhibit  a  plateau,  with  the  location  of  the  hump  related 
to  the  filter  cutoff  scale.  Thus,  the  beginning  of  the 
plateau  is  an  indirect  measure  of  the  largest  scales. 
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Fig.  4.  Longitudinal  velocity  structure  functions  as  a  function  of 
separation  distance  for  eight  individual  flight  segments:  1-6  Aug 
1999,  9.7  km;  2-23  Aug  2001,  9.3  km;  3-6  Jun  2000,  11.4  km 
(Wales);  4-17  Aug  2001, 10.6  km;  5-10  May  2000, 12.3  km  (Wales); 
6-6  Jun  2000,  9.9  km  (Wales);  7-26  Aug  1998,  6.6  km;  8-20  Feb 
1999,  10.9  km  (Japan).  Dashed  line  shows  r2/3;  dotted  line  shows 
r2/5.  Gray  lines  and  symbols  represent  data  at  scales  above  aircraft 
inertia  range. 

The  average  structure  functions  for  the  individual 
flight  segment  were  ensemble-averaged  for  all  129  levels 
and  for  various  subensembles,  weighted  by  the  number 
of  data  points  for  each  flight-segment  average. 

4.  Results 

a.  Structure  functions  for  individual  segments 

Longitudinal  structure  functions  Z)LL  for  eight  indi¬ 
vidual  level  flight  segments  are  shown  in  Fig.  4.  Note  that 
above  1  km,  the  curves  are  shown  in  gray  to  denote  the 
region  where  aircraft  trajectories  are  likely  to  influence 
scaling  exponents.  These  cases  were  chosen  to  span  the 
entire  range  of  turbulence  strength  in  the  3D  scaling 
range,  approximately  evenly  spaced.  The  strong  turbu¬ 
lence  cases  (curves  1  through  3)  all  have  extensive  3D 
scaling  ranges,  with  scaling  exponents  close  to  2h  (dashed 
line).  In  contrast,  the  weaker  cases  (curves  4  through  8) 
have  smaller  scaling  ranges  with  exponents  closer  to  0.4 
(dotted  line).  Only  curve  1  displays  a  terminal  plateau. 


Fig.  5.  Temperature  structure  functions  as  a  function  of  sepa¬ 
ration  distance  for  eight  individual  flight  segments.  Curves  1, 2,  and 
3,  see  Fig.  4;  9-26  Aug  1999,  12.5  km;  10-11  May  2000,  12.3  km 
(Wales);  11-18  May  2000,  13.8  km  (Wales);  12-24  Aug  2001, 
8.7  km;  13-6  Aug  1999, 12.2  km.  Dashed  line  shows  r2/3;  dotted  line 
shows  r2/5.  Gray  lines  and  symbols  represent  data  at  scales  above 
aircraft  inertia  range. 

while  curves  2  through  4  show  a  plateau  level  at  in¬ 
termediate  scales  (100  to  1000  nr),  followed  by  a  sub¬ 
sequent  increase  in  values  at  larger  scales;  this  behavior 
indicates  a  distinct  separation  of  small  and  large  scales. 
The  cases  with  weaker  small-scale  turbulence,  curves  5 
through  8,  do  not  have  a  plateau  level  but  instead  display 
increasing  slope  above  the  shrinking  3D  scaling  range. 

Temperature  structure  functions  are  shown  in  Fig.  5 
for  eight  flight  levels,  with  curves  1, 2,  and  3  representing 
the  same  levels  as  those  in  Fig.  4.  A  striking  feature  is  the 
wide  range  of  slopes  in  the  3D  scaling  range.  In  addition, 
only  curves  1  and  3  show  plateaus,  similar  to  their  cor¬ 
responding  /J|  [  curves  in  Fig.  4.  The  rest  of  the  cases 
display  increasing  slopes  above  the  3D  scaling  range. 

A  notable  feature  seen  in  Figs.  4  and  5  is  the  incon¬ 
sistent  relation  between  structure  functions  at  small  and 
large  scales,  best  seen  in  curve  6  in  Fig.  4  and  curve  10  in 
Fig.  5.  These  both  feature  weak  small-scale  turbulence, 
but  large-scale  turbulence  among  the  strongest  of  all 
cases.  (See  section  5a  for  more  discussion). 
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Fig.  6.  Vertical  velocity  structure  functions  as  a  function  of  sepa¬ 
ration  distance  for  eight  individual  flight  segments.  Curves  1-8,  see 
Fig.  4;  14-12  Sep  2002,  10.4  km.  Dashed  line  shows  r273;  dotted  line 
shows  r215.  Gray  lines  and  symbols  represent  data  at  scales  above 
aircraft  inertia  range. 


Figure  6  includes  vertical  velocity  structure  functions 
for  eight  individual  segments,  with  all  curves  except  for 
curve  14  representing  the  same  flight  levels  as  those  in 
Fig.  4.  With  the  exception  of  curve  8,  all  feature  a  pla¬ 
teau.  The  plateaus  in  curves  1  through  3  are  terminal, 
while  curve  4  displays  an  intermediate  plateau.  In  curves 
5, 14,  and  7.  the  plateaus  begin  at  kilometer  scales,  with 
no  plateaus  in  the  corresponding  Z)LL  curve.  Fifty-one  of 
the  129  levels  display  some  type  of  plateau. 

Note  that  for  curves  1  through  4,  the  plateaus  begin  at 
scales  smaller  than  the  plateaus  in  the  corresponding 
D ll  curves  (Fig.  4),  indicating  that  limiting  scales  for  the 
vertical  fluctuations  are  smaller  than  those  for  the  hor¬ 
izontal  fluctuations.  This  suggests  that  the  limiting  scale 
for  the  vertical  fluctuations  is  the  vertical  height  of  the 
CAT  layer,  an  idea  supported  by  the  fact  that  the  ma¬ 
jority  of  the  plateaus  seen  in  Dww  begin  at  scales  from 
50  m  to  1  km,  echoing  the  range  of  UTLS  CAT  layer 
heights. 

Plateaus  are  generally  associated  with  the  strongest 
turbulence,  including  the  20  levels  with  the  strongest  small- 
scale  turbulence.  It  is  not  clear  whether  this  connection 


Fig.  7.  Distributions  of  exponents  nx  for  velocity  and  tempera¬ 
ture  structure  function;  Dxx « r"x,  where  X  =  L,  T,  W,  or  N. 
Vertical  dashed  line  is  %. 


is  rooted  in  physically  relevant  dynamics  or  is  simply  an 
indication  that  plateaus  for  weaker  turbulence  are  masked 
by  uncorrelated  vertical  motions.  In  the  latter  case,  the 
presence  or  lack  of  a  plateau  is  simply  a  consequence  of 
the  relative  strengths  of  the  small-  and  large-scale  ver¬ 
tical  motions. 

b.  Structure  function  scaling  exponents  and 
structure  parameters 

Distributions  of  the  3D  scaling  range  structure  func¬ 
tion  scaling  exponents  (e.g.,  nL  where  DLL°c /■"*•)  are 
shown  in  Fig.  7.  All  four  distributions  have  peak  values 
and  mean  values  between  0.6  and  0.7,  with  mean  values 
of  0.61  for  nL,  0.69  for  nN  and  nT,  and  0.62  for  nw.  The  n, 
distribution  has  the  narrowest  distribution,  with  90%  of 
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the  levels  having  an  exponent  between  0.4  and  0.75.  The 
distributions  for  nN  and  nT  are  more  symmetric  and 
broader  than  that  for  nL,  while  the  nw  distribution  is 
broader  still  with  a  secondary  peak  near  0.88.  Although 
79  of  the  129  cases  (61%)  have  values  of  nL  close  to  % 
(0.55  to  0.8),  only  21  of  those  have  values  of  nN,  nw ,  and 
nT  all  within  that  range  as  well. 

The  range  of  values  in  structure  function  exponents 
for  /J|  |  matches  well  with  results  of  stratified-turbulence 
DNS  from  Werne  and  Fritts  (2000);  they  report  that 
the  exponents  for  DI  (  and  />IT  concentrate  between  2/s 
and  %,  with  variations  associated  with  location  within 
the  shear  layer  and  the  temporal  development  stage 
of  the  layer.  This  suggests  that  the  range  of  structure 
function  scaling  exponents  seen  in  the  aircraft  data 
may  be  a  result  of  similar  differences  in  the  location 
and  the  time  at  which  the  CAT  layers  were  sampled  by 
the  aircraft.  However,  further  analysis  of  measured  small- 
scale  structure  functions  and  direct  numerical  simulation 
(DNS)  is  needed  to  verify  this.  Although  the  2/s  value  is 
a  theoretical  result  for  n T  in  a  buoyancy  range  (Bolgiano 
1959),  the  corresponding  value  of  1  'A  for  nL  was  not 
seen  in  the  DNS  of  Werne  and  Fritts  (2000),  nor  is  it 
evident  in  the  experimental  data  reported  here.  As  such, 
it  would  be  erroneous  to  assume  that  the  lower  values 
are  associated  with  buoyancy  subrange  scaling  as  de¬ 
scribed  by  Bolgiano  (1959). 

Even  though  there  is  a  wide  range  of  power-law  expo¬ 
nents  for  the  individual  segments,  the  ensemble-averaged 
structure  functions,  for  the  overall  ensemble  and  the 
subensembles,  do  feature  exponents  close  to  z/3  for  all 
velocities  and  temperature  (section  4c).  This  is  qualita¬ 
tively  similar  to  observations  by  Duck  and  White  way 
(2005),  who  found  that  spectra  for  short,  strong  turbu¬ 
lence  segments  behaved  as  k  2  while  averages  of  several 
segments  fit  well  with  k  5/3.  However,  the  range  of 
scaling  exponents  seen  here  is  likely  not  a  short-segment 
effect,  since  no  correlation  is  seen  between  segment 
length  and  scaling  exponent  for  individual  segments.  In 
fact,  the  shortest  level  flight  segment  of  45.9  km,  (which 
is  still  over  500  times  longer  than  the  end  of  the  3D 
scaling  range)  features  nL  of  0.663,  while  for  the  longest 
flight  segment,  386  km,  nL  =  0.54.  The  reason  that  the 
ensemble  averages  have  scaling  exponents  near  %  is  that 
the  distributions  of  scaling  exponents  peak  close  to  that 
value,  so  that  most  of  the  levels  that  make  up  the  en¬ 
semble  behave  in  that  manner. 

Figure  8a  shows  the  distributions  of  structure  param¬ 
eters  C/2  and  Ct,  from  the  r2'3  force-fits  of  /.)r  [  and  D  \-\ 
in  the  3D  scaling  range  (section  4a),  using  a  logi0  bin  size 
of  0.55.  These  follow  a  lognormal  distribution  (lines  in  the 
figure),  a  behavior  predicted  by  Kolmogorov  (1941)  for 
dissipation  and  also  shown  by  FS05.  The  upper  scale  in 


lo9ioe 


lo9io(Cx) 


Fig.  8.  (a)  Distributions  of  C|  (m4/3  s~2)  and  Ct  (K2  m~2/3)  for 
3D  scale  range,  using  a  log1()  bin  size  of  0.55;  scale  on  upper  axis  is 
dissipation  values  corresponding  to  C2  curve  based  on  Eq.  (4). 
Asterisks,  C2;  crosses,  Cf- ;  dashed  line,  lognormal  fit  for  C2;  dotted- 
dashed  line,  lognormal  fit  for  Ct-  (b)  Distributions  of  Cjv/C2 
and  Cw/Ct  for  3D  scale  range  using  a  log10  bin  size  of  0.12. 
Asterisks  =  Cat/C£;  crosses  =  C|£/C£;  vertical  line  marks  iso¬ 
tropic  value  of  4/3. 

Fig.  8a  indicates  the  dissipation  values  corresponding  to 
the  C2  distribution  based  on  Eq.  (4).  The  means  and 
standard  deviations  of  the  log  distributions  are  as  follows: 
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=  0.719; 

(c2l) 

=  5.51  X 
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s-2; 

(6) 
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(J\  — 
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1(T4  m2  s 

-3. 

1 

(7) 

(l°gioCr> 

=  -4.29; 

CT,  ' 

loglOC7- 

=  0.749; 

(Ct) 

=  2.77  X 

nr4  k2  nr2/3. 

(8) 

The  distributions  of  the  structure  parameters  for  the 
normal  and  vertical  velocities  are  reported  in  Fig.  8b  in 
the  form  of  ratios  against  C'f  and  are  best  matched  by 
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a  lognormal  distribution.  The  peak  of  the  log  distribu¬ 
tion  for  CnICl  corresponds  to  a  value  of  1.2,  close  to  the 
isotropic  value  of  4/3.  However,  only  44  of  the  129  levels 
display  a  value  within  ±20%  of  the  isotropic  value.  The 
log  distribution  for  Cw/Cf  exhibits  greater  deviation 
from  isotropic  values,  centered  on  a  peak  corresponding 
to  0.40,  and  with  only  3  of  the  129  levels  showing  a  value 
within  ±20%  of  the  isotropic  value. 

The  significant  variation  of  the  lateral  structure  pa¬ 
rameters  from  the  isotropic  value  of  4/3  may  be  an  ex¬ 
planation  for  the  wide  range  of  scaling  exponents  seen  in 
Fig.  7  (i.e.,  greater  deviation  from  anisotropy  may  lead 
to  scaling  exponents  farther  from  2A).  However,  no  cor¬ 
relation  was  observed  between  Cw/C2  or  CnICl  and  the 
scaling  exponents  for  the  individual  levels. 

Because  of  the  wide  range  of  scaling  exponents  shown 
in  Fig.  7,  the  force-fit  values  of  structure  parameters 
should  not  be  confused  with  the  classic  structure  pa¬ 
rameter  defined  for  true  r2B  scaling  (i.e.,  the  value  of  the 
structure  function  at  r  =  1  m).  Rather,  they  represent  a 
useful  metric  for  comparing  relative  turbulence  strengths 
in  the  3D  scaling  range. 

c.  Ensemble-averaged  structure  functions  and 

parameters 

Ensemble-averaged  structure  functions  for  velocity 
and  temperature  as  a  function  of  separation  distance  are 
shown  in  Fig.  9,  and  the  scaling  exponents  and  structure 
parameters  are  included  in  column  1  of  Table  2.  The  3D 
scaling  range  for  the  velocity  extends  from  2  to  70  m, 
with  exponents  close  to  ZA  for  DLL  and  DNN  but  larger 
for  Dww.  The  scaling  range  for  temperature  is  between 
100  and  1000  m,  with  an  exponent  of  0.7.  Consistent  with 
statistics  of  individual  flight  levels  (Fig.  8),  the  normal 
and  vertical  velocity  structure  parameters  reflect  aniso¬ 
tropic  behavior,  with  Cf,!Cf  =  0.88  and  C'fvICj  =  0.6.  At 
larger  scales,  greater  than  about  8  km,  the  power-law 
exponents  of  DLL  and  /)  rr  range  are  between  0.75  and 
0.85,  although  these  values  should  be  viewed  with  caution 
because  of  the  small  size  of  the  linear  range  as  well  as  the 
likely  effect  of  aircraft  trajectories  at  these  scales  (see 
section  2).  At  intermediate  scales,  Z)NN  and  /Jr  [  show 
a  decreasing  slope  and  DT T  an  increasing  slope,  whereas 
Dw\\i  diverges  sharply  from  DLL,  the  latter  indicating 
a  transition  away  from  three-dimensional  behavior. 

Although  the  ensemble-averaged  structure  functions 
for  all  129  levels  are  useful,  they  reveal  only  part  of  the 
picture  because  they  are  weighted  heavily  toward  the 
strongest  turbulence  levels.  To  better  study  behavior  for 
all  turbulence  levels,  ensemble  averages  were  found  for 
three  subensembles,  with  results  shown  in  Fig.  10  and 
Table  2.  The  three  subensembles  for  velocity  structure 
functions  are  as  follows: 


Fig.  9.  Composite  structure  functions  as  a  function  of  separation 
distance  for  all  levels.  Solid  line:  £>LL  (m2  s-2);  dashed  line:  Z)NN 
(m2  s~2);  solid  line  with  plus  symbols:  DWw  (m2  s~2);  dotted- 
dashed  line:  (K2).  Thin  dashed  lines:  r2/3. 

•  Weak  turbulence:  [log10C|  —  (log10C^)]  <  -*1Ao_log  ci 

[see  Eq.  (6)], 

•  Moderate  turbulence: 

-1/2CTloSlfl  -  [loSl0 Cl  ~  (lQgl0 CL>]  <  1/2£7log10C£’ 

.  Strong  turbulence:  [log10C£  -  (log10C£)]  >  . 

Similar  subensembles  were  defined  for  Djt  ensembles 
based  on  values  of  log10(Cr). 

For  both  Z)LL  and  /JTT,  the  behavior  at  large  scales  is 
nearly  the  same  for  all  ensembles,  clearly  illustrating  the 
lack  of  correspondence  between  small-  and  large-scale 
values  mentioned  in  section  4a.  This  suggests  a  “decou¬ 
pling”  of  behavior  in  these  two  scale  ranges.  As  such,  the 
intermediate  range  scalings  are  mainly  functions  of  the 
strength  of  the  small-scale  turbulence,  with  higher 
slopes  for  weaker  turbulence.  Following  the  trend  seen 
in  the  overall  ensemble  averages,  the  3D  scaling  range 
exponents  are  close  to  2A  for  the  horizontal  velocities 
and  temperature,  the  exception  being  the  weak  turbu¬ 
lence  ensemble  with  values  for  nL,  nN,  and  nw  closer  to 
0.5.  The  ratio  CwICf  is  approximately  0.6  for  all  en¬ 
sembles,  while  CnICl  varies  between  0.82  and  1.0  with 
no  apparent  trend. 

5.  Discussion 

The  structure  function  data  presented  in  section  4 
revealed  two  important  features  that  will  be  explored  in 
more  detail  in  this  section — the  decoupling  of  the  small 
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Table  2.  Structure  parameters  and  isotropy  ratios  for  ensemble-averaged  structure  functions. 


log10C£  range 

All 

>-2.52 

-3.24  to  -2.52 

<-3.24 

nL 

0.66 

0.67 

0.71 

0.46 

nN 

0.70 

0.71 

0.73 

0.52 

nw 

0.76 

0.82 

0.6 

0.46 

Cl  (m4/3  s“2) 

6.22  X  10~3 

1.85  X  HT2 

1.44  X  10-3 

3.17  X  HT4 

s  (m2  s  3) 

1.73  X  10~4 

9.02  X  10-4 

1.93  X  10-5 

1.99  X  10~6 

cltcl 

0.88 

0.87 

1.00 

0.82 

cltcl 

0.60 

0.61 

0.64 

0.62 

logiuC£  range 

ALL 

>-3.92 

-4.67  to  -3.92 

<-4.67 

Cl  (K2  nL2'3) 

3.25  X  10~4 

9.26  X  10-4 

6.54  X  10'5 

1.10  X  10~5 

rif 

0.70 

0.78 

0.67 

0.67 

and  large  scales,  and  the  deviation  of  structure  param¬ 
eter  from  isotropic  values  .  The  focus  will  be  on  framing 
these  small-scale  behaviors  in  the  context  of  3D  strati¬ 
fied  shear  turbulence.  Such  turbulence  develops  when 
the  local  Richardson  number  falls  below  a  critical  value 
because  of  high  shear  or  weak  stratification  (Hazel  1972); 
the  Richardson  number  is  given  by  Ri  =  N2IS2,  where 
S  is  the  shear  rate  {=[(duLldz)2  +  ( duNldz)2]u 2}  and  N is 
the  buoyancy  frequency  {=[g/0(50/5z)]1/2}  where  9  is 
potential  temperature. 

Laboratory  experiments  (Miles  1961;  Thorpe  1973a; 
DaSilva  et  al.  1996)  have  shown  that  the  transition  to 
turbulence  is  via  the  development  of  2D  Kelvin-Helmholtz 
(KH)  billows  that  become  unstable  and  eventually  break¬ 
down  into  3D  turbulence.  Low  Reynolds  direct  numerical 
simulations  have  provided  additional  insight  into  the 
evolution  of  stratified  shear  layer  turbulence;  these  in¬ 
clude  the  work  of  Smyth  and  Mourn  (2000),  Smyth  et  al. 
(2001),  Joseph  et  al.  (2004),  Palmer  et  al.  (1994,  1996), 
and  Werne  and  Fritts  (1999,  2000). 

The  developing  KH  billow  grows  to  a  height  larger 
than  the  initial  shear-layer  thickness.  Shear  layers  with 
strong  initial  stability  (i.e.,  larger  Ri)  result  in  thinner 
turbulence  layers  relative  to  the  initial  shear  layer  height, 
with  vertical  velocity  levels  and  scales  limited.  The  KH- 
to-turbulence  transition  process  features  several  phases: 
2D  billow  formation  and  rollup,  initial  transverse  in¬ 
stability,  transition  to  turbulence,  and  breakdown  of  the 
billow  structure  leading  to  a  mainly  parallel  flow  layer 
with  decaying  turbulence  and  internal  gravity  waves. 
Scaling  and  statistical  characteristics  vary  among  these 
phases  and  among  different  vertical  locations  within  the 
layer  (e.g.,  Smyth  and  Mourn  2000;  Smyth  et  al.  2001; 
Joseph  et  al.  2004).  Unfortunately,  these  effects  are 
difficult  to  capture  with  aircraft  measurements,  since 
there  is  no  way  to  know  the  “age”  of  the  layer  or  the 
vertical  location  within  the  layer. 

Kelvin-Helmholtz  instabilities  are  not  the  only  possi¬ 
ble  source  of  the  small-scale  turbulence  described  here, 
but  they  will  be  the  focus  of  the  discussion  in  subsequent 


sections.  For  other  possible  sources,  see  Whiteway  et  al. 
(2003a,  2004)  and  Duck  and  Whiteway  (2005),  who  in¬ 
vestigated  gravity  wave  features  of  several  of  the  flights 
reported  here,  and  Riley  and  Lindborg  (2008),  who 


Fig.  10.  Composite  structure  functions  as  a  function  of  separa¬ 
tion  distance  for  subensembles  of  levels.  Solid  line,  all  levels  (same 
as  Fig.  9);  dashed  line,  weak  turbulence — all  levels  with  log(Ci)  or 
log(C£)  <  Via  below  mean;  solid  line  with  plus  symbols,  moderate 
turbulence — all  levels  with  log(C£)  or  log(Cr)  within  Via  of  mean; 
dotted-dashed  line,  strong  turbulence — all  levels  with  log(C|)  or 
log(Cr)  >  Via  above  mean.  Thin  dashed  lines  show  rm\  thin  dotted 
line  shows  r215. 


April  2010 


WROBLEWSKI  ET  AL. 


1167 


proposed  a  model  for  stratified  shear  turbulence  based 
on  the  breakup  of  large-scale  horizontal  motions. 

a.  Large-  and  small-scale  decoupling 

The  notion  of  decoupling  of  small-  and  large-scale 
behavior  refers  to  the  fact  that  the  strength  of  turbulence 
at  large  scales  is  not  necessarily  a  good  predictor  of  the 
strength  at  small  scales.  It  is  most  evident  in  the  individual 
level  structure  functions  in  Figs.  4-6  and  the  ensemble 
averages  in  Figs.  9  and  10.  It  can  also  be  seen  in  the  variety 
of  scaling  behavior  in  the  intermediate  scales  (1  to  10  km) 
that  may  not  necessarily  represent  any  physically  relevant 
mechanism  but  is  simply  a  “bridge”  between  the  two 
distinct  scale  ranges,  adjusting  to  their  relative  strengths. 
Given  this,  extrapolation  of  large-scale  behavior  is  a 
questionable  approach  for  inferring  small-scale  turbu¬ 
lence,  as  proposed  in  FS05.  The  most  striking  evidence 
of  decoupling  is  the  plateau  levels  observed  in  DWw; 
clearly,  vertical  velocity  is  the  best  indicator  of  transition 
between  3D  small-scale  turbulence  and  large-scale  mo¬ 
tions  dominated  by  horizontal  velocities. 

The  source  of  energy  for  the  small-scale  turbulence 
seen  in  the  aircraft  measurements  is  at  the  scale  of  the 
most  unstable  wavelength  of  the  initial  shear  instability, 
on  the  order  of  several  kilometers.  Larger-scale  motions 
that  may  be  present  would  not  necessarily  correlate  with 
the  stratified  shear  turbulence  at  smaller  scales,  leading 
to  the  decoupling  seen  in  the  structure  functions.  The 
lack  of  variability  of  the  large-scale  behavior,  best  seen 
in  the  ensemble-averaged  Z)L l  and  DTT  curves  in  Fig.  10, 
suggests  that  the  mechanism  responsible  for  these  larger- 
scale  motions  may  be  more  ubiquitous  and  uniform  than 
those  associated  with  the  thin  shear  layers,  which  appear 
to  be  more  intermittent  and  variable  based  on  larger 
spread  in  the  sub-kilometer  structure  functions.  In  es¬ 
sence,  the  large-scale  motions  develop  independently 
from  sub-kilometer-scale  CAT  layers.  The  exception  to 
this  would  be  situations  for  which  the  source  of  initial 
shear  is  connected  to  the  large-scale  motions;  for  exam¬ 
ple,  Riley  and  Lindborg  (2008)  conjecture  that  energy 
injected  at  large  horizontal  scales  could  break  down  into 
smaller  structures  in  the  presence  of  strong  stratification 
resulting  in  vertical  shearing,  development  of  3D  turbu¬ 
lence  layers,  and  a  downward  cascade  of  energy. 

For  the  majority  of  the  129  cases  presented  here,  the  jet 
stream  is  the  likely  source  of  the  initial  shear,  since  flight 
planning  protocol  focused  on  locating  and  targeting  strong 
winds  and  shear  associated  with  jet  streams.  The  maxi¬ 
mum  jet  stream  speeds,  based  on  atmospheric  soundings 
at  nearby  locations  (http://weather.uwyo.edu/upperair/ 
sounding.html),  were  over  50  kt  for  106  of  the  levels  and 
over  90  kt  for  98  levels.  As  a  way  to  test  the  importance 
of  the  jet  stream  as  a  source  of  shear  in  the  aircraft 
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Fig.  11.  The  quantities  C2  (m4/3  s-2)  (triangles)  and  Ct( K2  nT2/3) 
(circles)  as  a  function  of  vertical  distance  from  jet  centerline  z  zjet. 
Values  are  averages  over  0.65-m  bins  in  Z  Zjet. 

measurements,  the  levels  were  binned  based  on  abso¬ 
lute  distance  between  the  measured  CAT  level  and  the 
center  of  the  jet  found  from  the  sounding  profiles.  As 
seen  in  Fig.  11,  bin-averaged  values  of  C\  and  Cl  (with 
bin  sizes  of  0.65  km)  are  highest  at  vertical  locations 
between  2  and  4  km  from  the  jet  center,  decreasing  on 
either  side,  corresponding  to  high-shear  regions  along 
the  wings  of  the  jet. 

The  value  of  the  structure  functions  at  scales  greater 
than  about  5  km  are  likely  contaminated  by  the  trajectory- 
anisotropy  effects  discussed  in  section  2  (LST04;  Lovejoy 
et  al.  2009;  Lilley  et  al.  2008).  Although  the  analysis  of 
such  effects  is  not  the  focus  of  this  paper,  it  is  worth 
considering  whether  they  are  the  source  of  the  observed 
decoupling  between  small  and  large  scales;  that  is,  is 
the  decoupling  simply  an  artifact  of  aircraft  measure¬ 
ments  rather  than  a  physical  phenomenon  rooted  in  the 
fluid  dynamics  of  the  atmosphere?  Several  points  argue 
against  that  explanation,  all  based  on  the  fact  that  any 
artifacts  in  the  measurements  due  to  the  trajectory  ef¬ 
fects  should  appear  at  scales  greater  than  about  1  km, 
the  end  of  the  range  for  which  aircraft  inertia  smoothes 
out  trajectory  effects.  Figure  10  shows  that  the  dif¬ 
ferences  in  turbulence  strength  that  lead  to  the  idea  of 
decoupling  occur  at  small  scales,  with  the  large-scale 
structure  functions  being  nearly  the  same.  Even  if  the 
large-scale  behavior  is  primarily  an  artifact,  it  cannot 
explain  the  difference  in  turbulence  levels  at  small  scales. 
Along  the  same  line,  Lovejoy  et  al.  (2009)  showed  that 
an  initial  “break”  (change  in  scaling  exponent)  in  the 
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spectra  would  occur  at  the  end  of  the  inertial  range  of 
the  aircraft,  but  changes  in  slopes  of  the  structure 
functions  reported  here  occur  at  scales  as  low  as  tens  of 
meters  (see  curves  7  and  8  on  Fig.  4),  well  within  the 
aircraft  inertial  range.  In  addition,  Lovejoy  et  al.  (2009) 
showed  that  temperature  and  humidity  spectra  did  not 
display  the  same  breaks  as  horizontal  velocity,  yet  the 
decoupling  effect  and  the  changes  in  slope  are  seen  here 
in  temperature  as  well  (seen  most  clearly  in  Fig.  10). 
Finally,  the  most  obvious  evidence  of  the  decoupling  is 
the  plateau  features  of  the  vertical  velocity  structure 
function  that  signal  a  clear  separation  of  scales;  for  all 
but  2  of  the  51  cases  that  display  such  plateaus,  these 
features  begin  well  within  the  aircraft  inertial  region,  as 
determined  by  the  Dz  data,  and  therefore  cannot  be 
explained  as  a  trajectory-based  artifacts. 

b.  Anisotropy 1 * * 

Despite  the  deviation  of  C%  and  Cw  from  isotropic 
values  of  4/3C£  in  the  3D  scaling  range,  the  degree  of 
anisotropy  is  relatively  small  in  the  sense  that  ratios  of 
Cl/Cl  and  Cw/Cf  are  of  0(1).  As  such,  estimates 
of  dissipation  found  from  Eq.  (4)  are  probably  reason¬ 
able  in  an  order-of-magnitude  sense,  though  further  work 
is  needed  to  quantify  the  extent  to  which  this  is  true.  The 
degree  of  anisotropy  is  larger  in  Cw  than  in  C]^  because  of 
the  effect  of  buoyancy,  which  suppresses  vertical  mo¬ 
tions  but  not  to  the  extent  that  the  flow  becomes  two- 
dimensional.  Gargett  (1988)  relates  these  two  situations 
to  the  degree  of  stratification  and  would  classify  the 
weakly  anisotropic  condition  seen  here  as  “isotropic.” 

The  level  of  anisotropy  in  Cw  is  correlated  to  the 
strength  of  the  small-scale  turbulence,  with  greater  an¬ 
isotropy  for  weaker  turbulence  cases.  This  is  seen  in  the 
Cw/d  distributions  in  Fig.  12a  for  the  weak  and  strong 
cases  defined  in  section  4c,  with  the  distribution  peaking 
at  0.28  for  the  weak  ensemble  and  0.47  for  the  strong 
ensemble.  Figure  12b  further  illustrates  the  anisotropy 
trend,  showing  bin-averaged  values  of  C^/Cl  increasing 
with  increasing  C£. 

DNS  and  laboratory  experiments  of  stratified  shear 
turbulence  have  shown  that  Ri  increases  from  its  initial 
shear  layer  value  as  the  layer  develops,  leveling  off  at 
values  ranging  from  0.25  to  0.4  during  the  latter  phases 
of  layer  life  when  turbulence  levels  decay  rapidly  (Smyth 
and  Mourn  2000;  Thorpe  1973b).  These  stronger  levels 
of  stratification  give  rise  to  internal  gravity  waves,  and 
if  the  weaker  turbulence  cases  reported  here  represent 


1  Note  that  anisotropy  here  refers  to  deviation  of  structure  pa¬ 

rameters  from  the  isotropic  values  and  is  not  an  indication  of  an¬ 

isotropic  scaling  as  described  by  Lovejoy  et  al.  (2007). 


Fig.  12.  (a)  Distributions  of  ratio  C\y/C£  for  3D  scale  range  using 
a  log10  bin  size  of  0.15.  Asterisks  indicate  strong  turbulence  (see 
Fig.  10  caption);  crosses,  weak  turbulence.  Dashed  vertical  line 
represents  isotropic  value  of  4/3.  (b)  The  ratio  Cj£/C£  as  a  func¬ 
tion  of  C£.  Values  are  averages  over  bins  of  logi0(C£),  with  a  bin 
size  of  0.6. 

latter-stage  shear  layers,  then  their  higher  levels  of  anisot¬ 
ropy  may  also  be  a  result  of  the  stronger  stratification.  In 
support  of  this,  structure  function  analysis  of  stratified- 
shear  DNS  revealed  C(v/C£  =  0.43  in  the  midlayer  for 
the  stage  of  evolution  characterized  by  decaying  turbu¬ 
lence  (Werne  and  Fritts  2000).  This  value  is  close  to  those 
seen  in  Fig.  12b  for  the  weakest  turbulence  cases. 

6.  Concluding  remarks 

High-resolution  measurements  obtained  from  a  unique 
aircraft  platform  have  shed  light  on  turbulence  behavior 
in  the  sub-kilometer  range  for  which  published  data  is 
scarce.  At  scales  less  than  100  m,  structure  functions 
exhibited  power-law  exponents  that  ranged  from  0.4  to 
0.75  for  individual  flight  segments,  but  values  close  to  % 
for  the  overall  ensemble  and  for  various  subensembles. 
Weakly  anisotropic  behavior  was  seen,  with  structure 
parameters  Cn  and  Cw  on  the  same  order  as  C£  but 
deviating  from  the  expected  isotropic  value  of  4/3 C£. 
Anisotropy  in  Cw  was  stronger  on  average  for  weaker 
turbulence  levels,  consistent  with  DNS  findings.  This  3D 
subrange  was  decoupled  from  behavior  at  scales  greater 
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than  10  km,  which  exhibited  less  variation  in  the  structure 
function  values,  suggesting  that  mechanisms  associated 
with  small-scale  stratified  shear  turbulence  are  more  in¬ 
termittent  and  variable  than  those  associated  with  large- 
scale  behavior.  The  observed  decoupling  calls  into  question 
the  use  of  mesoscale  data  to  infer  microscale  behavior.  The 
T>ww  curves  displayed  plateau  levels  characteristic  of  a 
limiting  vertical  scale,  a  feature  generally  associated 
with  stronger  turbulence  levels  in  the  3D  scaling  range. 
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